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Abstract: Strong correlation of photons, particularly in the single-photon 
regime, has recently been exploited for various applications in quantum 
information processing. Existing correlation measurements, however, do 
not fully characterize multi-photon correlation in a relevant context and may 
pose limitations in practical situations. We propose a conceptually rigorous, 
but easy-to-implement, criterion for detecting correlated multi-photon 
emission out of a quantum optical system, drawn from the context of 
wavefunction collapse. We illustrate the robustness of our approach against 
experimental limitations by considering an anharmonic optical system. 
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1. Introduction 

Strong correlation of photons at the few quanta level can make possible a variety of nonlinear 
optical devices useful for quantum information processing, such as single-photon transistors or 
switching devices mill and the generation of single photons on demand |[3]|4l. Furthermore, the 
photon-photon correlation mediated by the emitters can also be employed to simulate quantum 
many-body systems in a controllable way. For example, the effective on-site repulsion between 
photons can be exploited to study quantum phase transitions such as the Mott-superfluid tran- 
sition El |6] 111 m and the fermionization of bosons fOl. These applications are closely related 
to a specific correlation effect, namely a single-photon blockade effect ||2l [3] [TO] [TTI — Once 
a system is excited by one photon, the abosorption of next photons is blocked, e.g., due to the 
anharmonic energy level structure of the system. 

Recently, interest in the correlation effect has also been extended to multi-photon level in 
the context of multi-photon gateway, where a random (Poissonian) stream of photons can be 
converted into a bunch of temporally correlated n photons. In particular, Kubanek et al. demon- 
strated the operation of two-photon gateway to some extent using an optical cavity QED system 
lfT2l . In view of all these efforts, it seems very crucial to have a theoretical framework that can 
appropriately characterize multi-photon correlations lfT3l[T4l[T5l[T6l . e.g., «-photon blockade 
effect, and desirably that can be efficiently tested in experiment. 

Conventionally, correlation of photons is measured by the nth-order coherence functions 
introduced by Glauber [iT], g(")(0) = (fl^"fl" )/{a'a)'\ where a (a') is the annihilation (cre- 
ation) operator of an optical field. However, it is noted in the recent experiments of cavity QED 
Il2]|3][l2 that §'^'(0) is not effective to resolve the correlated two photon emission due to a 
huge bunching at the atom-cavity bare resonance overshadowing the two-photon resonance. 
Furthermore, as n goes beyond two, g^"^ (0) or its simple variants [11211 contain more peaks at 
k = 1, •••,«— 1 photon resonances, irrelevant to genuine «-photon correlation, as to be shown 
below. 



Instead one may take n-photon excitation peaks in {a^"a'^) spectrum itself, e.g. in lfT2l[T3l[T4l 
[15] US), as a confirming evidence of «-photon correlations. Rigorously speaking, however, the 
multi-photon resonant excitation peaks spectroscopically identified only uncover the energy- 
level structure of the system. Whether each peak in the bare coincidence {a' "a") indicates 
relevant photon correlation must be checked very carefully. For example, the three-photon co- 
incidence tends to increase in the spectrum without any correlation if the system possesses 
two-photon correlation since an uncorrelated emission added to a correlated pair may register 
another three-photon coincidence. Therefore, we need to consider a stricter physical context for 
characterizing correlated emission of «-photons. Although there have been several studies on 
higher-order photon statistics in view of nonclassicality lfT8l[T9l . none of them carries a clear 
interpretation as multiphoton correlation. 

Here we propose a conceptually rigorous, but easy-to-implement, criterion for measuring 
correlated multi-photon emissions out of a quantum optical system. The criterion is derived 
by considering wavefunction collapse related to sequential photo-detections ||20| and has the 
following merits, (i) It only detects highly-correlated «-photon emissions with no classical ana- 
logue, and (ii) can be tested by usual photon coincidence detections (no conditional measure- 
ment) in a detector-efficiency insensitive form, and therefore, experimentally favorable. We 
also introduce a quantitative measure of «-photon correlation based on this criterion, which is 
quite robust in addressing multi-photon correlations against experimental imperfections. This is 
a practical merit of our approach made possible by the well-established correlation context. We 
illustrate the power of our method in an optical cavity QED system, where genuine n-photon 
correlated emissions can be efficiently verified in accordance with its anharmonic energy levels. 

2. The criterion 

In order to envision a generic, though not exhaustive, scenario where multiphoton correlations 
may arise, let us compare two systems, one with harmonic and the other with anharmonic 
level structure [Fig.l (a)]. When a harmonic system with level spacing haiQ is excited by an 
external driving on resonance (Oz, = (Bo), all energy levels are equally accessible. On the other 
hand, for an anharmonic system, if the external field is «-photon resonant with the «th level, 
other levels than «th would not be substantially addressed by the external field. As a result, 
the system could be excited to contain only n correlated quanta and further excitation would 
be prohibited — «-quanta (photon) blockade effect. We will apply a similar line of reasoning to 
emission, rather than excitation, process. Specifically, we construct a criterion to detect 'pure' 
n-photon correlated emission by incorporating two distinct features, (i) surge or rapid emission 
of photons up to n quanta and (ii) blockade beyond n, which can be applied to any quantum 
optical systems, not necessarily anharmonic ones. 

2.1. Photon surge 

Generally, the photo-detection rate is proportional to the intensity of the optical field un- 
der consideration, M °^ where the operators S'± correspond to the positive- and the 
negative-frequency part of the field. Let us assume that a quantum system can be described by 
a pure steady state \^)s for simplicity, but our argument applies equally well to mixed states. If 
it has emitted n — 1 quanta, the wavefunction is collapsed to |^[" = (f"^' , conditioned 
on these emissions. The detection rate for the succeeding nth photon is then given by 



n — 



{i'li'D 



(1) 



after the normalization of the conditional state \^[" ''). Specifically, if the emission out of 
the system is a bunch of highly correlated n-photons, the second photon will be emitted right 
after the first photon and the third photon after the second, and so on. This idea can be used to 
construct our criterion as follows. 

The "bare" rate for the first emission is simply given by the intensity, = {S_S^, which 
only characterizes the signal strength and has little to do with correlation. For «-photon correla- 
tion (n > 1), once a photon is emitted, however, the next emission will immediately follow, thus 
the conditional rate ^2 must be large enough. In particular, we require ^2 to be larger than !M\ , 

i.e. R2 1 = ^ = > 1, which is nothing but the bunching condition in the Glauber g^^^i 

function. Extending the requirement to next emissions sequentially, we derive a set of surge 
conditions 

which must be satisfied for each k = 2, . . . ,n. 

2.2. Photon blockade 

However, the satisfaction of Eq. (|2]i for all k ~ 2, . . . ,n is not sufficient to ensure «-photon 
correlation, and importantly, one must also look at the next occurrences carefully. After the 
detection of « photons, the succeeding emissions must be suppressed, which can be expressed 
as 

R-ti-i<l, ik = n + l,---). (3) 

The fulfillment of all the surge and the blockade conditions in Eqs. dU and (O respectively 
constitutes our criterion for «-photon correlated emission. Note that the condition Q coincides 
with the special case of the higher-order antibunching criteria introduced in ITS l for the non- 
classicality of photon statistics, rather than the correlation effect . 

In our criterion, it is crucial to use the conditional rates ^i^, rather than the bare rates ), 
as the former takes into account the correlation between adjacent emissions in a stronger sense. 
However, the resulting criterion does not require any conditional measurements. Instead, the 
quantities Rk,k-i in Eqs. ^ and (O simply involve various photon-coincidence rates and we 
particularly note that the numerator and the denominator are in the same order of the field 
strength. It is thus given in an experimentally desirable form, that is, insensitive to the quantum 
efficiency of photodetectors. 

2.3. Measure of multi-photon correlation 

The conditions in Eqs. (|2]i and ([3j may be used to define a quantitative measure ^„ of n-photon 
correlation as 

n Mr 

^„ = nmax{Ri.n-l,0} max{R,/_i - 1,0}, (4) 

k=2 k=,i+l 

where A'tr is a truncated excitation number to be taken appropriate to a given situation. ^„ 
quantifies the strength of the n-photon correlation by measuring the deviation of Rk.k- 1 from 
unity in the surge and the blockade conditions of Eqs. (|2]i and Q, respectively, and returns a 
nonzero value only when all those conditions are satisfied. To experimentally obtain ^„ for 
a given system, one first measures the bare A;-photon coincidence rates {S'^S'^l) for all k = 
1, . . . jA^tr. Then, each conditional rate Ra,a:-i defined by Eq. (|2j is evaluated and plugged in to 
Eq. to determine the value of 



2.4. Remarks 



(a) Conventionally, multi-photon correlations have been discussed in terms of the Glauber co- 
herence functions 

(„) _ (i_ (^1 (X2) ■■■i- {Xn)i+{Xn) ■ ■ ■ i+{x2)S+ (xj )) 
^ " {i^{x,)i+{x,)){i^{x2)S+{x2))---{i-{Xn)S+{Xn)) 

(jc,: a general space-time point) iflTl . The context of correlation in ^'"^ is, however, rather limited 
and we particularly note that g'"' compares the n-photon coincidence rate (numerator) only 
with the single-photon counting rates (denominator). Large (small) value of g'"' characterizes 
a bunching (antibunching) effect with no strict n-photon correlation that can emerge even in 
a classical scattering system, e.g. g^"^ = n\ for a thermal light (Hanbury-Brown-Twiss effect 
Il2ni22l ). Another example of g'^^ ^ 1 with no rigorous two-photon correlation will be shown 
below in Sec|3] 

(b) It is, therefore, interesting to ask whether our criteria of «-photon correlation can be 
fulfilled by a classical source. It turns out that, as mentioned in Sec. 12.21 the blockade condition 
in Eq. ^ is related to the nonclassicality of light fields ifTsl [191 . Let us consider the single- 
mode case in which the field amplitude ((f_) may be replaced by the annihilation (creation) 
operator a (a^). Then, as a special case of Ref. ifTSl . one can show that, for a classical source 
represented by a positive-definite Glauber-P function, P{a) > 0, a Cauchy-Schwarz inequality 
follows as 

The violation of the above inequality, which is nothing but the blockade condition Rk.k- 1 < 1 
in Eq. (|3]l, is thus a clear signature of nonclassicality. So our criteria of multiphoton correlation 
can be fulfiled only by nonclassical sources. We emphasize that, in the so called multiphoton 
antibunching criteria in ifTSl [191 , the focus was made on how to reveal nonclassicality of the 
field by a mathematical approach based on the positive Glauber-f function, thus lacking a clear 
interpretation as multiphoton correlation. 



3. Application: Cavity QED system 

3.1. Model 

To illustrate our criterion, we consider a cavity QED system — one of the well known anhar- 
monic systems that can be implemented in various experimental platforms ll2l[3l[4l [T4l . A qubit 
(two-state atom, quantum dot, etc.) is coupled to a single mode field driven by a classical field. 
For simplicity we investigate the on-resonance case, (Oa = COc = COo, where coa is the qubit tran- 
sition frequency and (Oc the cavity resonance frequency. The qubit-cavity system at coupling 
strength g is then described by the Hamiltonian 

H = h(0o ^fl^fl + ^CT;^ + ihg{a^a- - ao+), (6) 

where a± and Ut are the Pauli pseudospin operators. The composite system has the ground state 
|0,^) with the energy £0 = and the polaritonic excited states = (|n,g) ± |« — I,e)) 

with E„,± = iiTicoo±Tigy/n {n ~ 1,2...) [See Fig.l (b)]. Therefore, when the system is driven 
by an external field at frequency (Ol, «-photon resonant absorption may occur f23\ at 



nhcOi = £';i,± = nhcOQ ± hg\/n. 



(7) 




Fig. 1. (a) Energy-level diagram for (i) harmonic and (ii) anharmonic system, (b) Energy 
level structure for cavity QED system, (c) multiphoton coincidence rates (a'"fl") as a func- 
tion of 8/g for 2k/ g = y/g = 0.01 with S /k = 0.1. The dotted vertical lines represent the 
locations of the multiphoton resonances, 8 = ^gls/n throughout Figs. 1-3. 



In practical situations, the qubit and the cavity field may interact with Markovian envi- 
ronments, which causes dissipation and decoherence to the system. The global evolution is 
then governed by the master equation p = ^[///,p] + 7(<7_pa+ — ^a+a-p — jpcT+CT-) + 
K'(2flpfl^ — a^ap — pa^a), where 7 {2k) is the qubit (cavity) decay rate and the interaction 
Hamiltonian Hj = hd [a^a + 5(7,) + ihg{a'0- — aa+) + ihS{a^ — a), with the driving strength 
(o and the detuning 5 = (Oq — (Ol- 

By measuring the cavity transmission as the driving frequency (Ol scanned, one may identify 
the energy-level structure of the cavity QED system. In Fig. [Tfc), we plot the bare «-photon 
coincidence rate, as a function of the normalized detuning 5/^. In the weak-excitation 

limit, these rates are related to the n-excitation probability P„ as {a^"a") « n\P„. We see that 
more resonant peaks are spectroscopically observed at 5 = ±g/y/n as the order n is increased 
in a very strong-coupling regime, y/g = 2K/g = 0.01. It is important to note that not all the 
peaks in nth order coincidence rate are relevant to «-photon correlation (e.g., the peak at 5 = 
zLg in the two-photon coincidence), so the bare coincidence rates may not be used as such to 
address genuine multiphoton correlation. To overcome this difficulty, for instance, one may 
try to classify those peaks with a prior knowledge on the excitation paths ITSl . However, in 
realistic situations, the peaks become less resolved as the coupling strength is reduced (not 
shown). More importantly, these resonant peaks give information only on the level structure of 
excitation and have a weak connection to correlated emissions. 

3.2. Correlation measures 

Instead, if one measures the Glauber coherence function g'^"^ of the output, the result may char- 
acterize the correlation of emitted photons to some extent, but not in a full rigorous sense. 
In particular, ^'"'(0) = {a''"a")/{a^a)" in Fig. |2a) shows a large bunching at zero detuning 
5 = 0, which has nothing to do with genuine «-photon correlation as we will clearly show 
below. Close inspection of photon statistics reveals that the system does exhibit some non- 
classical behavior at 5 = 0, e.g. the oscillation of conditional detection rate which peaks 
at even number of k, but it is not a rigorous «-photon correlation at any level n in view of 
our criterion. To get rid of this "cumbersome" resonance effect observed at 5 = that may 
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Fig. 2. (a)-(c) The conventional correlation functions g'"' (0) by Glauber and C^"' (0) by 
Kubanek etal. 1121 . (d)-(f) The quantitative measure ^„ together with conditional relative 
rates K/c^k-l- The truncation numbers used are (d) N^^ = 4, (e) and (f) = 5. In all plots, 
2K/g = Y/g = <^/K = 0.1. 



overwhelm the other resonance peaks, Kubanek et al. introduced the differential correlation 
function, C'^^(O) = {a''^a^) — (a^a)^, that measures the absolute occurrence of two-photon 
excitation with respect to the single-photon excitation 1121 . The context in this correlation 
function, however, is insufficient just like g^^\0) in general, although it was instrumental to 
identify the second resonant peak in flU . Furthermore, a generalization to «-photon level, 

for « > 3, becomes hardly effective in identifying the higher-order 
peaks by the broadening effect in the realistic regime [Fig. 2 (a)-(c)]. 

In contrast, our criterion not only detects correlated emission in a well-defined context, but 
also provides a practical tool to identify the multi-photon resonance structure of a system in 
realistic situations. In Figs.|2d),|2{e), and|2f), we plot the quantitative measure of Eq. ^ 
for n = 2,3, and 4 along with various rates Rk.k-i which are ingredients for constructing the 
corresponding y£'„. Remarkably, ^„ yields a positive value only in the spectral vicinity of the 
resonant peaks 5 = ztg/^/n. The "spurious" peak at 5 = disappears by our criterion, which 
rigorously confirms that this seeming "resonance" indeed does not represent pure «-photon 
correlated emission. 
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Fig. 3. Comparison between ^'■"'(0) (black dotted curve) and our measure ^„ (red solid 
curve) for (a) n = 2 and (b) n = 3, with the truncation numbers (a) A'ti- = 4 and (b) = 5, 
respectively. The driving intensity is rather high, S'/k = 1, with the coupling condition 
2K/g = Y/g = 0.1. 



3.3. Large driving field 

As we increase the pump strength to obtain more substantial signal, the coincidence spectrum 
usually becomes difficult to resolve due to the saturation of the system. Our method is, how- 
ever, still useful for moderately strong pumping owing to the rigorous context established in 
it. To demonstrate this merit, we have considered the case of a large driving field S Ik= 1 in 
Fig.|3] together with the realistic coupling yjg = 2fc/g = 0.1. Due to the intensity-dependent 
broadening effect, the Glauber function g'"^ (0) no longer shows noticeable marks of resonance, 
except for the peak at 5 = overwhelming the entire shape in the spectrum. On the other hand, 
our measure ^„ identifies a clear signature of multi-photon correlations under the same condi- 
tion. This capability would allow one to increase the pump strength to some extent, and thereby 
easing the difficulty of having to measure higher-order coincidence than n [i.e., blockade con- 
ditions in Eq. (|3]l] to identify «-photon correlation in our method. Furthermore, we have also 
checked that other possible broadening effects, e.g. atomic motion in the cavity, do not de- 
grade the capability of our criterion for characterizing multiphoton correlations. We attribute 
this robustness against experimental imperfections to the rigorous context established with the 
measure 



4. Conclusion 



In conclusion, we have devised an easy-to-implement criterion for detecting correlated multi- 
photon emission, imposing surge and blockade requirements in photoemission processes. A 
quantitative measure ./#„ has been derived from the correlation context between successive pho- 
ton emissions in the framework of wavefunction collapse. Our criterion applies to any quantum 
optical systems, including the ones with anharmonic structure (cavity QED systems, multi-level 
atoms, etc.). 

We have illustrated our method can efficiently detect multi-photon correlations at the reso- 
nant peaks of the cavity QED system in contrast to the existing correlation functions. Note that 
the anharmonic spectrum which scales as ^/n is a clear signature of quantum nature of light 
field 12411251 . and it thus has been of considerable interest for long but experimentally verified 
only recently ifTSl [T4l [TSll . In the optical cavity-QED system ll3ll4l [T2l[T3]| . it becomes harder 
to directly observe this anharmonicity in higher-orders due to less strong coupling than in the 
microwave circuit QED system, but our method remarkably makes it possible to clearly pick 
up the Y^-dependence despite experimental limitations. We anticipate that our conceptually 
rigorous approach can also be useful in addressing correlation effects in other quantum systems 
beyond optics. 
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